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ABSTRACT 



We present the results of a Cosmic Microwave Background (CMB) anisotropy 
ground-based millimetric experiment sensitive to fluctuations on angular scales 
of ~ 2°. Four independent bands centered at 3.3, 2.1, 1.3 and 1.1 mm collected 
~ 550 hours of observation during the Summer of 1994. The instrument 
was located on the island of Tenerife at an altitude of 2400 m. The low 
water-vapor content and the atmospheric stability of the site, combined 
with new techniques to subtract atmospheric noise, result in the reduction 
of atmospheric contamination in the lowest frequency channel to a level of 
~ 1.5 times the instrument noise. Detailed estimations of Galactic foreground 
contamination show that this contribution is neghgible at |6| > 12°. Two 
different multipole bands {i = 53llil ^ind 331^3) are analyzed showing that 
our technique to subtract the atmospheric contribution is more effective in 
the multipole band at ^ = 53. A hkelihood analysis of these data reveals the 
presence of a common signal between the channels at 3.3, 2.1 and 1.3 mm 
corresponding to a band power estimate of + l)Q/(27r) = 2.0lJ:^ " lO'^ 
and ^Je{e+ l)Ci/{2n) = 4.ll|^ • 10"^ for the £ = 53 and 33 multipole bands 
respectively. Calibration uncertainty has been treated as a systematic effect. 
The level of fluctuations in the £ = 53 band is in good agreement with our 
preliminary analysis presented in Piccirillo et al. 1997, with measurements by 
other experiments working at similar angular scales, and with the predictions of 
standard Cold Dark Matter (CDM) models. 



Subject headings: Cosmology: cosmic microwave background - Observations 
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1. INTRODUCTION 



In recent years the COBE DMR ([Bennett et al. 1996|) , FIRS ( |Ganga et ai 1993|) 



and Tenerife ([Hancock et al. 1994|) experiments have successfully detected anisotropy 



at large angular scales yielding precise estimations of the overall normalization of the 
power spectrum ( [Hancock et al. 1997a[ ) and its shape at low multipole moments £ ^ 30. 



At smaller angular scales {I ^ 70) a wealth of experiments have reported detections 
of anisotropy (for an up-to-date list of experimental results see Tegmark's web site: 
http: / /www.sns.ias.edu/^max/ cmb/experiments.htm][ ). These detections clearly suggest a 



increase in the power spectrum at I ~ 200 ( [Netterheld et al. 1997 , R£i-iicock et al. 1997b ) 



and its decline at ^ ~ 400 ([Scott et al. 1996] ). The main aim of our experiment (for 



preliminary results see piccirillo et al. 1997[) is to fill the gap in £-space remaining between 



the two £-space regions described above. Thus, with our observing strategy we are probing 
the multipole bands at i = 33ti3 ^ — 531^3, in between the Tenerife experiment at 
£ = 20 ± 8 ( [Hancock et al. 1997a[ ) and the ACME South Pole at i = 6811 ([Gundersen et 
al. 1995[) . 



The choice of the observing frequencies also distinguishes this experiment from other 
ground-based experiments. It is known that the combined contributions from Galactic 
foregrounds and discrete radio sources reaches its minimum at about 55 GHz. However, 
the presence of a strong O2 absorption line in the range 50-70 GHz, and the increasing 
atmospheric emission at millimetric wavelengths hamper observations at these frequencies 
forces ground-based experiments to observe at cm wavelengths where contamination from 
synchrotron and free-free emission are of concern. Thus, the millimetric range has been 
traditionally exploited by satellite and balloon experiments, and only very few ground-based 
experiments have observed at these frequencies (Python (piatt et al. 1997 ), SuZIE ( [Church 
et al. I997D and [Andreani et al. 1991[) . 
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In section ^ we briefly describe the instrument and concentrate on the measurement 
technique. In section ^ we give the details of the observations. Atmospheric effects in our 
data are described in section ^. The details of the atmospheric noise reduction technique 
and the resulting final data sets are presented in section ^. In section |] we conclude that at 
high |6| the only expected signal is the CMB. Section |^ describes the statistical data, and 
conclusions are presented in section ^. 



2. INSTRUMENTAL SET-UP AND MEASUREMENT TECHNIQUE 



The instrument has been described in detail in previous work (Piccirillo 1991| ; [Piccirillo 



fc Calisse 1993|) . The optics forms an off-axis Gregorian telescope with a parabolic 



primary mirror (45 cm diameter) and a hyperbolic secondary mirror (28 cm diameter). 
The detector is a four-channel photometer equipped with ^He bolometers working at 
0.33K. The bands peak at 3.3, 2.1, 1.3 and 1.1 mm wavelengths (channels 1, 2, 3 and 4 
respectively). The beam response for all channels can be well approximated by a Gaussian 
with FWHM=2.°03 ± 0.°09 and no significant sidelobes are found. 

In addition to the optics, the response of the instrument to a point-like source depends 
on the observing strategy and demodulation of the data ( [White fc Srednicki 1995| ). The 



observing strategy consists in daily drift scans at constant declination achieved by fixing 
the telescope in azimuth (0) and elevation {6). Additionally the beam moves on the sky as 
the primary mirror wobbles sinusoidally while the secondary is fixed. The right ascension 
(a) and declination (6) at which the antenna is pointing at time t are given in a good 
approximation by: 6(t) = 6o and a{t) = + (3q/ cos(5o) x sin(27r/^t + e), where (ao,5o) 
is the initial position of the antenna, (3q = 2.6° is the zero-to-peak chopping amplitude at 
a reference frequency = 4 Hz and e is an irrelevant phase constant. Each bolometer's 
output is sampled at fs = 80 Hz coherently with the mirror movement so to have 20 
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samples per mirror oscillation. The signal is demodulated in software by evaluating the 
amplitude of the first (4 Hz) and second (8 Hz) harmonic of the reference frequency fw 
Each demodulation (IF and 2F for the first and second harmonic respectively) can be 
divided in two components: the in-phase component containing mostly the sky-signal 
coherent with the reference motion plus the instrument noise and the out-phase containing 
mostly the instrument noise plus other sources of systematic noise. Such division requires 
a careful choice of phase which was obtained from observation of the Moon transiting the 
instrument beam. Then the mapping function at a sky location of coordinates {a, 6) when 
the antenna is pointing towards {ao,6o) is computed as ( [White fc Srednicki 1995|) : 

Nn /-i/Ca/u,) 
MnF{a,S;ao,6o) = — - / cos(n 27r/^t + () exp 

VSttct^ J-l/{2U) 

for n = 1,2 indicating respectively the IF and 2F demodulation, a = FWHM/ VS In 2, ( 
is a phase constant, A{t) is the angular distance between the point of coordinates {a, 6) and 
the center of the beam at time t and is a normalization constant for the n demodulation 
obtained by requiring an output of 1 K when the input is an extended source of 1 K filling 
the positive lobe. 

The response of the IF and 2F demodulations resembles the usual 2-beam and 3-beam 
responses to the transit of a point-like source through the beam. The IF demodulated 
data are well fitted by a 2-beam response with asymmetric Gaussians with cTq = l.°03 in 
the RA direction and = 0.°86 in the declination direction, and with a beam throw of 
Pq = 2.°38. The 2F demodulated data are fitted by a 3-beam response with aa = l.°56, 
(Ts = 0.°86 for the positive lobe, and cTq, = 0.°89, as = 0.°86 for the negative lobes and 
a beam throw of (3q = 2.°40. Additionally, these fits must be multiplied by the factors 
'RiF = 0.362 and = 0.593 to yield the normalizations to the actual response functions 
for IF and 2F respectively. These approximations greatly simplify the statistical analyses 
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as discussed in section |^. 

The response of the instrument to different angular scales (represented in terms of the 
corresponding multipole moment t) is completely specified by the window function. For a 
constant declination observation the window function for the product of two temperatures 
separated by an angle if) is computed according to [White fc Srednicki 1995| to give: 



{I - 2r)A$ 



(2) 



where Bi refers to the beam profile: Bi{a) = exp[—i{i + l)cr/2], J„ is the n-th order Bessel 
function of the first kind, jo is the zeroth-order spherical Bessel function and A$ is the 
bin size in radians on the sky. In figure ^ we show the window functions for ip = and 
A$ = 4 min in RA, corresponding to the ^-ranges iip = 33^13 and = ^•^-il- The 
central value corresponds to the band power average and the upper and lower limits give 
the i ranges where We{0)n has amplitudes larger than e^^^^ times its peak. 



3. CALIBRATION AND OBSERVATIONS 

Laboratory calibrations were performed by placing blackbody radiators at different 
temperatures in front of the optical window. By means of an off-axis mirror the optics 
within the cryostat is redirected towards a vessel containing eccosorb and divided in four 
sections. Two sections were filled with liquid Nitrogen alternately placed between the other 
two sections which were filled with liquid Oxygen. A measurement of the pressure at which 
these two liquids evaporate gives a precise measurement of the temperature. Then the 
container is rotated at 2 revolutions per second so the detectors see two black bodies of 
known temperatures. The calibration factors thus obtained have uncertainties ranging from 
1.5% (Channel 1) up to 7.6% (Channel 4). Since we are concerned with the possibility of 
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systematic effects in tlie calibration process we quote a laboratory calibration uncertainty 
based on our observations of the Moon (see below). 

Our observations were conducted during June and July 1994, collecting about 550 
hours of data, at Observatorio del Teide at Tenerife (Spain). The site, at an altitude of 
2400 m and latitude N 28°29, is well known among the solar community for its good seeing 
and the large percentage of clear days, the latter due to the fact that the inversion layer lies 
below the observatory for about 75 % of the time. This is also the location of the Tenerife 



experiment ([Hancock et al. 1994|) which has already demonstrated the potential of the site 
for hosting observations in the cm range due to the stability of the atmosphere and its 
excellent transparency ( Davies et al. 19^ ). These features makes of this site a promising 



place for mm observations. 

The observations concentrated at declination 5 = 40° by pointing the antenna towards 
the North meridian and at an elevation angle oi h = 78°. 7. This declination has been 
extensively measured from this site at larger angular scales(~ 5°) and lower frequencies 
(10, 15 and 33 GHz) with reported detections of structures in the CMB by the Jodrell 
Bank-IAC experiments ( [Hancock et al. 1994| ; |Gutierrez et al. 1997|) . 



The measurements started after the Sun was well below the horizon in order to 
avoid solar contamination. The angular distance between Moon and beam was always 
greater than 23°. 5 for the 5 = 40° observations. In addition the Moon was observed for 
astronomical calibration. The values obtained from the analysis of the Moon transits in 
the IF demodulation agree within ~ 20 % of the laboratory calibrations for all channels 
(see figure |^). The bulk of the calibration uncertainty (6 %, 26 %, 22 % and 10 % for 
channels 1, 2, 3 and 4 respectively) we believe is due to the lunar models used ( |Bennett~e^ 
al. 1992a[ , [Hagfors 1970| ). By fitting the observed Moon transit to the expected transits at 



both demodulations, and modeling the Moon as an extended uniform disk and the beam 
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as Gaussian shaped, we obtain estimates for the beam width, cr, and chopping amphtude, 
/3o, given in table ^ Since the Moon transit is better defined in the IF data we adopt the 
estimates for a and (3q from the IF demodulation. The values for /5o from the fits to each 
channel are consistent with each other, yielding an average value of (3o = 2°. 90 ± 0°.03. 
As expected, the beam widths seen by each channel show a monotonic decrease with 
frequency so that channel 1 has the widest beam and channel 4 the narrowest beam. Even 
so, given the sizes of the error bars assigned to each width, all of them are consistent 
with a unique beam width given by the weighted average: a = 0°.86 ± 0°.04. In any 
case, for calibration purposes we adopted for each channel the best fit value obtained for 
that channel. It is worth mentioning the agreement between the FWHM estimated for 
channel 1 in the lab {FWHM = 2°. 40 ± 0°.10) with that from the fit to the Moon transit 
{FWHM = 2°.21 ± 0°.24 from the IF Moon transit and FWHM = 2°.4 ± 0°.6 from the 2F 
Moon transit). The discrepancy between the values of (3o measured in the lab and from the 
Moon transit can only be attributed to the antenna assembling, and throughout this work 
we have assumed the value of (3q as obtained from the Moon transit analysis. 

4. ATMOSPHERIC EFFECTS ON THE DATA 

The biggest source of noise for CMB ground observations at millimetric wavelengths 
comes from atmospheric emissions dominated by fiuctuations of O2 and water vapor 
contents. Given the stability of the pressure during the campaign, the most important 
parameter turned out to be the precipitable water vapor. The meteorological conditions 
during the campaign are summarized in figure ^ The pressure and temperature are 
measured at the observatory four times a day. The relative humidity is also measured at 
the observatory every few minutes. The precipitable water vapor is obtained from the 
measurements of the balloons launched twice a day by the Spanish Meteorological Institute 
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from sea level. In our case, the precipitable water vapor is obtained by integrating the 
balloon measurements from the Observatory level up to 12 — 15 Km. 

4.1. ATMOSPHERIC NOISE LEVELS VS INSTRUMENT NOISE LEVELS 

Atmospheric effects in our data become evident as a big increase with respect to 
the expected instrument noise levels. In the ideal case (i.e. the instrument noise is 
white) we should observe flat power spectra if our data were entirely due to instrumental 
noise. However our data are a combination of instrument noise, atmospheric noise and 
astronomical signal, the latter to be ignored given its weakness in each night of observation. 
At low frequencies neither the instrument noise nor the atmospheric noise exhibit fiat 
spectra. Bolometers show a characteristic 1// spectrum while atmospheric fluctuations 
have a more complicated spectrum and in general follow a l/f^ spectrum with k variable. 
On the other hand, both spectra flatten at high frequencies, becoming indistinguishable 
from each other. At high frequencies the power spectra are flat, and assigning to the 
instrument noise the value corresponding to such noise floor is a huge overestimation. 

In theory, this problem could be solved by using the out-phase component produced 
during the demodulation because it is expected to contain only instrumental noise. In 
practice, even in the out-phase component there are still residual amounts of atmospheric 
noise because it is impossible to flnd a constant demodulation phase which completely 
sets to zero the sky-signal in the out-phase component. In flgure ^ we show the noise 
spectra in thermodynamic temperature for all channels and both demodulations of the in- 
and out-phase components for a typical night of observation. Channel 1 exhibits power 
spectra almost flat in the out-phase spectra for IF and 2F demodulations, indicating 
that most of the noise in channel 1 at high frequencies is due to instrument noise. The 
out-phase components for the rest of the channels still contains considerable amounts of 
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atmospheric noise as indicated by the similar shapes of the spectra in the in- and out-phase 
components. This interpretation is strongly supported by the fact that the spectral shape at 
low frequencies, and the noise floor at high frequencies, change from day to day. The higher 
values for the instrumental noise as obtained from the analysis of the out-phase components 
of the the IF data indicate the greater ability of the 2F demodulation in removing linear 
gradients caused by atmospheric emission . The final upper limits in thermodynamic units 
assigned to the instrument noise are 2.8, 0.8, 1.5 and 1.0 mK s^/^ for channels 1 to 4 as 
obtained from the analysis of the out-phase components of the 2F demodulation. In figure ^ 
we show the distribution of the rms values, calculated over 30 minute intervals, with 10 
second bins for all channels and both demodulations. At 10 seconds the contribution from 
instrument noise to the rms would be less than ImK for all channels, so most of the noise 
must be atmospheric. The noise level and the width of its distribution are larger in the 
IF demodulation. Additionally, the minimum values obtained in the 2F data are smaller; 
another indication of the enhanced ability of the 2F demodulation to reduce atmospheric 
noise. 



4.2. ATMOSPHERIC CORRELATIONS IN THE DATA 

The second relevant effect due to atmospheric contamination concerns correlation 
between different channels. This is evident from the highly correlated time-variable signals 
seen by all channels for the same night of observation. This point will be used subsequently 
on as the cornerstone for our analysis to reduce atmospheric noise in the data. In table ^ 
we present the mean correlation between channel i {i = 1, 2, 3) and channel 4 for the whole 
campaign and for the data selected to build the final data sets. As expected this correlation 
increases when increasing the frequency of channel i and it is higher in the IF data than in 
the 2F data. 
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Another effect is the auto-correlation introduced by atmospheric noise. Primary 
evidence of this effect is that the rms of the data do not follow the 1/\/N law when the 
data are binned. We have computed the auto-correlation function for the data within 
the same night of observation as seen by each channel at both demodulations. This was 
done by using the power spectra of the in-phase components and the Wiener- Khinchin 
relations (see e.g. [Barkat 19911 ). As noticed when computing the instrument noise limits, 
the in-phase component power spectra show a large departure from white noise, implying 
non-zero auto-correlations at lags different from zero. The average auto-correlation function 
for the whole campaign and for all channels and both demodulations is shown in figure ^. 
For a given lag we observe that the auto-correlation function is always larger in the high 
frequency channels. The coherence time is always larger for the IF data than for the 2F 
data. This reflects the higher efficiency with which a double-switching technique is able 
to discard the effects of linear temporal drifts due to atmospheric emission. From figure |^ 
we conclude that it is necessary to use bins large enough such that correlation between 
adjacent bins is reduced (see section This will ensure that the standard deviations are 
properly computed and assigned as error bars when moving to larger bin sizes. 



5. DATA PROCESSING 

The demodulation process produces a data point every 0.25 s per each IF and 2F 
demodulation and per channel. During the demodulation 3.3% of the data was rejected due 
to problems of synchronism between the mirror movement and the data acquisition system. 
The bulk of the data rejection was performed during the subsequent phases of binning and 
editing comprising the data processing. After demodulation a binning is performed to bring 
the data from 0.25 s to 10 s. This is a trade-off between the need to reduce correlation 
between adjacent bins (auto-correlation < 65% for all channels at IF demodulation and 
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< 35% for all channels at 2F demodulation) and to have a large enough number of points 
for the linear fits and further binning processes involved in the cleaning technique. During 
the 0.25 s to 10 s binning an iterative 3cr filter is applied three times to discard glitches due 
to malfunctions of the data acquisition system. Bins of 10 s built with 15 or less points at 
0.25 s are also discarded. In total, the amount of data rejected in this binning is about 8% 
for all channels and both demodulations. 

As explained in the section ^, the extremely large temporal correlation between all 
channels and the amplitudes and behavior of the noise levels with the bin size strongly 
support the atmospheric origin of the bulk of the noise in our data. Hence, the primary goal 
of any data reduction concerns the assessment and subtraction/reduction of this unwanted 
source of noise. 

5.1. DESCRIPTION OF THE CLEANING TECHNIQUE 

The approach adopted in this work consists in exploiting the high correlation between 
channels as an indicator of the atmospheric emission, so that by using a channel as a monitor 
we can clean the rest of the channels. Since channel 4 (1.1 mm band) is the most sensitive to 
atmospheric emission, we adopt it as the monitor. The method assumes that at each channel 
i we have a superposition of astronomical signal attenuated by fi due to the atmospheric 
transparency, plus the contribution from the atmospheric emission. When expressing all 
quantities in antenna temperature we have at channel i: ATant,! = fi^T^j^^rf^^ + AT^^ j. 
We assume that the atmospheric contributions are perfectly correlated between different 
channels so: AT^^vTi — ^i^'^ANTi- This assumption is strongly justified given the high 
correlation between data at different channels taken during the same night of observation 
(see table The amplitude of the signals are too high to be attributed to CMB signal or 
other astronomical signal (see section ^ giving additional support to the above assumption. 
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To recover the astronomical signal at channel i in thermodynamic temperature (AT"''*''°) 
and referred to the top of the atmosphere we solve the linear equation: 

ATant, = ^AT;-*- + {ATanta - ^—ATr'n x a, (3) 

Q C4 Pi4 

where ATANT,i and ATant,4: are the recorded antenna temperatures at channels i and 4 
respectively; fi and the atmospheric attenuation factors; Cj and C4 are the Rayleigh- Jeans 
to thermodynamic conversion factors:ATAiVT,i = -ATi {a =1.29, 1.66, 3.66 ,4.82 for 
channels 1 to 4 respectively) and ^,4 is the fraction of the astronomical signals seen at 
channels i and 4. The second term within the parenthesis accounts for common structure in 
channels i and 4, though with different amplitudes. Based on our estimations of non-CMB 
contaminants, at high Galactic latitude( high we expect no astronomical signal other 
than the CMB itself (see section H) so pj4 = 1. At low |6| and in the Galactic plane we 
expect different signals in shape and in amplitude as seen by different channels. This 
motivates the introduction of pi4 7^ 1, whose value is obtained from our estimations of the 
Galactic diffuse emission. 

The atmospheric cleaning technique is applied to the lOs-binned data after subtraction 
of a constant offset from each channel. The rejected data typically correspond to sections 
taken during bad weather and/or warming of the cryostat. The data are divided into 
segments of 5 minutes. For each of these segments we compute the values of from a linear 
fit of channel i versus channel 4. Likewise, all points within the same 5 minute segment 
share the same fi and /4 as computed from the splined values of water vapor w, pressure 
P and temperature T from figure ^ used as input to the code by [Cernicharo 1985| , which 



computes the atmospheric opacities due to water vapor and oxygen using the US standard 
atmosphere model. 
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5.2. PRODUCING THE FINAL DATA SETS 

After cleaning, a new binning is performed to bring each processed night of observation 
from a 10 second to a 4 minute bin size so the beam is sampled with at least 3 points. Once 
again a 3cr filter is applied to discard possible glitches occurred during the cleaning process. 
To each point at 4 minutes we assign an error bar given by the standard deviation of the 
mean of all points at 10s within the bin at 4 minutes. A similar analysis to that performed 
in section ^2 allows us to obtain the mean auto-correlation function for the processed data 
in each channel and both demodulations. This analysis shows that for points separated by 4 
minutes the auto-correlation functions for all channels 1 to 3 and both demodulations range 
between 9% (Ch 2 IF) and 1% (Ch 1 2F). Therefore, our 4 minute bins can be considered 
uncorrelated as well as their error bars . In both demodulations and in all channels we 
observe residual baselines of very long periods. We proceed to remove these remnants by 
fitting linear combinations of sinusoidal functions after a re-edit of the data. The re-edit 
discards noisy sections which may affect the fitting process. The minimum period of the 
sinusoidal functions is chosen to be large enough so to remove signals corresponding to 
angular scales bigger than the ones to which the instrument is sensitive. Thus, for the IF 
data the minimum period is 90° in RA and a minimum period of 72° in RA for the 2F 
data. In figure ^ we display these various stages of the cleaning technique for a typical night 
as seen in all channels and both demodulations. In the last column of table ^ we give the 
mean amplitude of the baseline fits. Columns 1 and 2 show the percentage of total data 
used with respect to the original data at 0.25 s and the number of nights used to generate 
the final data sets. The percentage of data used is bigger in the 2F data as well as the 
number of used nights, with the exception of channel 3. The final data sets were obtained 
by stacking all individual baseline-cleaned nights where the rms did not exceed 0.65, 1.3 
and 2.5 mK for channels 1, 2 and 3 respectively in the IF demodulation and 0.4, 0.4 and 
2.0 mK for channels 1, 2 and 3 respectively in the 2F demodulation. Column 3 in table ^ 
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gives the mean rms for the surviving nights once residual basehnes have been removed. 
The stacking process consists of computing weighted averages, where to the ith bin in the 
final data set we assign a value and error bar given by: 



where the indices ij refer to bin i in night j, ^Tij and aij are the data point and standard 
deviation at bin i in night j, and rii is the number of nights used for this ith bin in the final 
data set. The final data sets in the regions before and after the Galactic Plane crossing are 
shown in figure |[ 



5.3. PERFORMANCE OF THE TECHNIQUE 

The efficiency of the atmospheric reduction process is best demonstrated by looking 
at the power spectra of the data before and after its application. In figure ^ we show the 
power spectra of both demodulations and for channels 1, 2 and 3 for a typical night before 
and after cleaning. This figure corresponds to the same data as in figure ^. The reduction 
in the noise level is evident from these plots, so that the corrected Ch 1 and Ch 2 in both 
demodulations approach the levels expected from instrument noise. We also notice the 
flattening, approaching the ideal behavior of white noise. For channel 3 there is also an 
overall decrement in the power spectra of both demodulations indicating that a substantial 
fraction of the atmospheric noise has been subtracted. However, the levels of the cleaned 
data for channel 3 still show residual atmospheric contamination. In table ^ we give the 
values attained by the noise spectrum in the in-phase components at different frequencies 
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and for both demodulations before and after cleaning for the same observing night as in 
figure ^. The corresponding values are also given for the out-phase component spectra. 
At low frequencies (i.e 0.001 Hz) where atmospheric effects are more evident the cleaned 
file exhibits lower values of the noise spectra than before cleaning; a good indication that 
most of the atmospheric noise has been removed. Further, this value approaches the values 
attained in the out-phase components containing small amounts of atmospheric noise due 
to the leakage during the demodulation. 

Another indication of the performance of the technique is given in figure where 
we show the distribution of the rms values measured in the 10 second binned data after 
applying our cleaning method. Once again the last bin exhibits a high value since it also 
contains the contributions from all following bins to avoid a huge spread in the plots. When 
we compare this figure with figure ^ we notice a huge decrement of the rms for all channels 
and both demodulations. 

Another check of the performance is the recovery of the Galactic Plane (GP) crossing. 
The introduction of a factor 7^ 1 to recover the GP explicitly assumes perfect correlation 
of this signal as seen by channel i and channel 4. This is a very good approximation for 
channels 2, 3 and 4 where the bulk of the Galactic emission is due to dust emission. The 
estimated GP crossings in channel 2 IF demodulation (2F demodulation) is 99% (96%) 
correlated with channel 4, while the estimation of the GP in channel 3 IF demodulation 
(2F demodulation) attains a 100% (100%) correlation with respect to the GP seen by 
channel 4. For channel 1 the contributions from dust and free-free emission are comparable, 
the latter slightly more important. Furthermore, the free-free template (1420 MHz map) 
and the dust template (240 pm DIRBE map) show a relative slight displacement in the 
position of the GP at 6 = 40°, thus lowering the correlation between channel 1 and channel 



4 down to 93% and 72% for the IF and 2F data respectively. In figure 11 we show the 
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predicted Galactic Plane crossings superimposed on our measurements for channels 1, 2 
and 3 and both IF and 2F demodulations. As discussed in section ^, we have considered 
the dust model in [Boulanger et al. 1996| , with the factors listed in table ^. We have also 



tried several different Galactic emission models in the literature and have checked that, 
while the absolute amplitudes in each channel depend strongly on the model used, their 
ratios are quite stable and so is the parameter pa as can be seen in table ^ The slight 
differences between the /9j4 factors for the IF and 2F demodulations are not relevant, except 
for channel 1. This is a direct consequence of the above mentioned displacements between 
the contributions generating the GP in channel 1. Using p[^ in equation (|^) instead of p^ 
results in an amplification/attenuation of the restored GP given by: 

p ^ fi/Cj - O^i ■ fi/{Ci ■ pa) 

fi/ci- ai- fi/{ci- p[^) 

The extreme values of Fj obtained when using the values obtained from the models 
described in table ^ are 0.6 and 1.2. Therefore we do not expect big changes in the 
amplitude of the restored GP due to the change of dust model used to estimate the values 
of the factors. The general agreement between the predictions and our measurements 
constitutes an important check on the performance of our system and method. 



6. GALACTIC FOREGROUNDS 

The Galactic contribution has been analyzed considering the synchrotron, free-free 
and dust emission. When studying the synchrotron emission it is convenient to distinguish 
between the expected levels of signal at low galactic latitudes (i.e at |6| < 12° or GP 
crossings) and at high galactic latitudes (i.e at |6| > 12° or outside the GP). 



As studied by p?egmark fc Efstathiou 1996| , outside the GP and at our angular scales 
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and frequencies, the contributions from the Galactic synchrotron and free-free emission are 
expected not to be a source of confusion with CMB anisotropy. We have also conducted 
a study of the expected levels of contamination by synchrotron and free-free emission in 
our data by convolving our instrumental response with a template of synchrotron and 
free-free processes. This template was obtained by extrapolating in frequency from the low 
frequency map at 1420 MHz ( [Reich fc Reich 19861 ), and modeling both processes with a 
single power-law: Tff^sync ■ The spectral index (3 is obtained by fitting this law to the 

low frequency maps at 408 MHz ([Haslam et al. 198^ and at 1420 MHz. The results of such 
an analysis predict that free-free plus synchrotron would contribute to the observed rms 
with values smaller than 1 fiK for both demodulations and in all channels. Additionally, 
assuming the extremely conservative approach that the signal seen at 33 GHz by the 
Tenerife experiment is entirely synchrotron emission, and extrapolating with (3 = 2.7, we 
obtain that for all channels the synchrotron emission at high |6| amounts for less than IfiK 
in both IF and 2F data. These two independent approaches clearly indicate the negligible 
effect of synchrotron and free- free emission on our data at high 

In the GP region the map at 1420 MHz is dominated by free-free emission due to the 
presence of the Cygnus X HII region and many other unresolved HII regions ( Davies 



al. 1996|). Accordingly, when generating the template for this region we assumed a spectral 



index of /? ~ 2.16 as obtained from equation (3) of [Bennett et al. 1992 b| . To generate 



the template in this region we also considered the possibility of still having non-negligible 
contributions of synchrotron emission in the 1420 MHz map as indicated by the spectral 
indexes required to reproduce the GP's seen by DMR at 31.5, 53 and 90 GHz. 

We have used the DIRBE map at 240 fim as template for the dust emission and 
extrapolated to our frequencies by using the fit to the FIRAS data obtained by [Boulangei| 
et al. 1996 resulting in a model with n = 2 and = 17.5K. We have checked the validity 
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of this model at low |6| by reproducing the expected dust contribution to the GP seen by 
COBE DMR at 90 GHz where the relevant section has been smoothed to 10° FWHM. 
Together with a brief description and reference to the dust model used, in table |^ we give 
the rms values expected from synchrotron plus free-free and dust emission for different dust 
models outside the GP > 12°). These rms values are completely negligible as compared 
to the observed rms values in our final data sets (see section [7p. This can also be seen by 
using the figure of AT^ust = 2.7 ± l.SfiK at 53 GHz and 10° angular resolution by p:logut 
al. 1996b|. 



Of much less concern is the contamination due to unresolved sources and from known 
point-like sources. The latter contribution was computed by convolving our instrument 
beam with a grid where we placed sources extracted from the Kiihr et al. (1981) and 
the Green Bank sky survey (|Gondon, Broderick fc Seielstad 1989|) complemented by the 
Michigan and Metsahovi monitoring programme. The weakest considered source presents 
a flux density at 5 GHz of 0.18 Jy. We extrapolate the flux density to our frequencies 
using the fit obtained by [Kiihr et al. 1981] where available. Fluxes of sources not present in 



Kiihr et al. 1981| , but for which we have measurements at three different frequencies, were 
fitted to a power law and extrapolated to our frequencies, while fiat spectra were assumed 
for those sources for which flux densities were available only at a single frequency. For all 
channels and both demodulations the expected rms in the section of our data |6| > 12° is 
much smaller than 1 fiK. Finally we refer to [Franceschini et al. 1989| to also exclude the 
contribution from randomly distributed sources given our observing frequencies and beam 
width. In view of all these figures we conclude that outside the GP the only expected 
astronomical signal should be CMB. 



7. STATISTICAL ANALYSIS 
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7.1. 



THE LIKELIHOOD FUNCTION 



The data selected for the statistical analysis have been chosen under two requirements. 
First the selected bins must correspond to regions of high \b\ where we are confident that 
Galactic emission is negligible. Thus we have identified two sections in our data: RAi 
for which b > +12° and RA2 for which b < —12°. The second requirement concerns 
the number of independent nights {rii in equations (H) and (|D) used to generate the final 
data sets. This causes the different RA ranges (columns 1 and 5 in table ^ over which 
the selected data for each channel and demodulation span. This happens because the 
different editions and cleaning technique are performed independently for each channel 
and demodulation. The minimum number of independent nights required is ni= 6, except 
for channels 1 and 2 in the 2F demodulation where the higher number of used nights 
allows to increase this threshold up to 9 and 8 points respectively. In all cases the number 
of points used for generating the final data bins is large enough to allow to use the 
Kolmogorov-Smirnov test to check that our bins in the final data sets are consistent with 
being drawn from Gaussian distributions. Then, the likelihood function for a set of data x 
is completely specified by giving the covariance matrix V and the vector of means x, which 
in our case is identically zero. Then: 



where |V| denotes the determinant of V. There are two independent contributions to the 
covariance matrix so that V = + V^i , where corresponds to the correlations 
between bins according to the model we are testing and shows dependence on the parameters 
to be estimated, and V/? is the data covariance matrix computed directly from the data. We 
compute the elements ij of Vt by assuming an intrinsic Gaussian auto-correlation function 
(GACF): Yintri^ij) = Co exp{—9fJ{2 9^)) where 9ij is the angular separation between bins 




(7) 
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i and j, 6c is the coherence angle, the angle of maximum sensitivity for each experimental 

1/2 

configuration: 6^ = 2.17° and 6c = 1.42° for the IF and 2F data respectively and Cq 
is the parameter to be estimated. We also have to consider the effects introduced by the 
finite width of the beam and the observing strategy. After convolving with the combination 
of Gaussian beams reproducing our instrument responses (see section we obtain that 
the elements of the matrices V^t,if and \^t,2F, for the IF and 2F demodulation 
respectively, are given by: 



Vt,if(^, j) = ^If X [2 ■ CM{6ij\ o-„, a^, as) - CM{6ij + /5; ^a-, (y&) - 

CM{6ij - P; cr„, cr„, as)] (8) 

^T,2F{i,j) = ^If>^ [CM{6ij;ac,crc,as) - CM{6ij + P;ac,(Ji,as) + 
-CM{6ij]ai,ai,as) — CM{6ij — (3;(Jc,(Ti,as) + 
^(Ca/(% + 2/5; ac, ai, as) + CM{6ij + 2/5; a„ au as))] (9) 

where Cm(^; a„ au as) = C^6l/^{61 + al + aj){6l + 2a|) exp[-^V(2(^c' + o^c + ^?))] 
and, as obtained in section ||, (cTq, a^, /5, Kii?) = (l.°03, 0.°86, 2.°38, 0.362) and 
(crc,ai,a5,/3,K2F) = (l.°56, 0.°89, 0.°86, 2.°40, 0.593) for the IF and 2F demodulation 
respectively. 

In a Bayesian interpretation with uniform prior, the likelihood (£) as a function of 

1/2 

the positive definite parameter Cq is directly proportional to the probability density 
function of Cg^^. Then the best estimation of C^^'^ is that value for which C is maximum 
while confidence levels [(Cd^^)i, {Cq^'^)2\ to a C% level have been computed by requiring 
C{{Cl'\) = Ci{C',%) and J^^?'/J' C^^^^^ = C/100. In 

what follows detections are given to a 68% CL, and upper limits to a 95% CL . A claim of 
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detection is made whenever the lower hmit of the confidence interval at 68% CL is not zero, 
otherwise we quote the upper limit at 95% CL. 



7.2. LIKELIHOOD ANALYSIS ON SINGLE CHANNELS 

We have applied the likelihood analysis to each channel and demodulation in the 
ranges: RAi, RA2 and RAi + RA2. The results of these analysis are given in table |^. In 
all cases the values obtained in all these RA ranges are consistent with the presence of a 
common signal. For channel 2 at IF demodulation the results between RAi and RA2 are 
still marginally consistent, but the result in the RAi section shows a strong dependence on 
the choice of individual nights to generate the final data set. Accordingly in what follows 
we will only consider the RA2 section for channel 2 at IF. 

For channels 1 and 2, in all valid ranges and in both demodulations, the signals detected 
are consistent between them and with values of ~ lOO/iK. These values are consistent with 
our previous results reported in [Piccirillo et al. 1997| , where we concluded that the slight 



excess of signal seen in channel 2 with respect to that obtained in channel 1 may indicate 
the presence of some residual levels of atmospheric noise at the same level as the expected 
CMB signal. 

The results obtained for channel 3 at both demodulations clearly indicate that it is still 

affected by important atmospheric residuals: CMB signal does not scale with frequency 

and our estimated signals due to diffuse Galactic contamination at channel 3 amount to 
1/2 

Co,Gai < 35 fiK, too large a difference to be caused by the uncertainties in the Galactic 
estimation procedure. 
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7.3. JOINT LIKELIHOOD ANALYSIS 

The joint analysis of all three channels allows us to estimate the most likely signal 
which is common to all of our frequencies, such as CMB anisotropy. To estimate the 
correlation between channels we have computed the cross-correlation function between 
sections of data at 10 s which overlap in the same nights in channels i and j. We make 
use of the generalized Wiener-Khinchin relations for stationary processes to obtain one 
cross-correlation curve per night. The average of these curves are plotted in figure O. We 



observe that the cross-correlation becomes negligible at scales smaller than our binning in 
the final data sets, being only significantly different from zero at zero-lags as indicated in 
columns 2 and 4 in table ^. This cross-correlation at zero-lag enhances the diagonal terms in 
the sub-matrices which take into account the correlation between different channels in the 
covariance matrix V/j. This effect has been analyzed and discussed in detail by [Gutierrez 



19971 , concluding that the net effect is an increase of the error bars as compared with the 
case where not such correlations are present. In table ^ we present the results from the 
joint analysis of any two channels and all three channels for each demodulation . These 
results have been obtained by using the whole data set except for those involving channel 2 
at IF for which only the RA2 section of the data set was used. 

The analysis on any combination of 2 channels indicates the presence of common signal, 

1/2 

which for the IF demodulation is Cq ~ ISO/xK and ~ 75/iK for the 2F demodulation. 
Although being completely consistent, these figures must be viewed with caution: the signal 
monotonically increases as we increase the frequencies of the channels being combined. 
This behavior again indicates the higher level of contamination in channel 3 for both 
demodulations. 

Finally, we have also considered the case of having a superposition of CMB signal plus 
a signal with a spectral behavior different from a black-body. In this way we obtain the 
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contamination due to a foreground component which is consistent with the data. We also 
assume a GACF with the same coherence angle as the CMB signal for the foreground signal 
in addition to a scaling of the signal with frequency. Thus the foreground signal in channel 
1 exhibits an intrinsic ACF Cci\{0) = Co,Fgdexp[-9'^ / while channel j (j = 2,3) 
shows Cctjie) = Ccl\{e) X X (cj/ci)2 . As in equation (D q (i = 1,2,3) is 

the Rayleigh- Jeans to thermodynamic conversion factor; Vi = 95.1 GHz, z/2 = 169.0 GHz 
and = 243.5 GHz the frequencies for channels 1, 2 and 3 respectively. By setting the 
value of n to the appropriate values we obtain the signal due to any of the foreground 
contaminants. The likelihood function now becomes a function of two parameters to be 
estimated: C = ^{Cl^Jmb^Cl^^pg^). Since, a priori, we do not have any information about 

1/2 

either of them we obtain the probability distribution function of Cq\,„^ by marginalizing 
with respect to C^pg^ and vice-versa. In table |10| we present the results of this analysis for 
different values of the spectral index n covering the ranges expected for dust, atmospheric. 



synchrotron and free-free emission. In figures ^ and |1J we show the contour plots for the 
spectral indexes n = 2.0, 0.0, -2.1 and -3.0, corresponding to dust, a simple atmospheric 
model, free-free and synchrotron emission in Rayleigh- Jeans approximation. The data are 
then converted into thermodynamic units in the plots. Small departures from these nominal 
values yield essentially the same results as seen in table 10. 

7.4. DISCUSSION OF THE RESULTS 

The joint analysis of channels 1, 2 and 3 at the 2F demodulation reveals the presence 
of a common sig nal with Cg^^ = 721^6 ^k. The conversion to band power estimates 



e.g. ISteinhardt 1994|) yields Ji{i + l)Q/(27r) = (2.0 ± 0.7) ■ 10"^ at £ = 53^1 



m 



good agreement with our previous results in [Piccirillo et al. 1997| (Cq^^ = 76^21 /uK 



for a GACF or Je{e + l)Ci/{27i) = (2.1 ± 0.6) ■ 10"^ in the band power estimate 
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notation). At similar angular scales, £ = 561^8, [Nettertield et al. 19971 quote a value of 

-5 



We have also tested the possibihty that part of the detected common signal is due to 
contamination by Galactic foregrounds or residual atmospheric contamination. The analysis 
confirms that dust can not be responsible of the detected signal (Co ^)^^,^ < 25 /iK), leaving 



a CMB signal of Cq^^ = 72^1 /iK which corresponds to + l)Q/(27r) = 2.0tE]:^ ■ 10"^ 
After considering calibration uncertainty as a systematic effect the figures above become 



^ 72+34 ^i{e+l)Ci/{2n) = 2.0tli ' 10"^ Less conclusive is our analysis 

when a free-free or synchrotron spectrum is assumed for the contaminant signal. Then we 
only have vague upper limits for the contaminants which were known in advance from both 
our estimations and results from the Tenerife experiment when observing at S = 40°. As 
indicated in section ^ the free-free contamination seen by the Tenerife experiment at 33 
GHz is about 4 /iK. The extrapolation in frequency renders this source of contamination 
to values < l/iK at both demodulations and in all channels. These upper limits are 
further reduced by the fact that our experiment probes higher values of ts than those 
probed by the Tenerife experiment, and at high |6| the Galactic power spectrum scales as 
d oc £-3( [legmark fc Etstathiou 1996| ; |Kogut et al. 1^^ ; [Gautier et al. 1992|) . 



The presence of atmospheric residuals is also tested by allowing for the presence of a 
signal with a spectral index in antenna temperature of n = 0., 0.2 and 0.4. The n = case 
corresponds to the approximate case in which the effective atmospheric temperature is the 
same in all our channels. The n = 0.2 and 0.4 cases allow an increase with frequency of the 
effective atmospheric temperature. For all the considered n values, the likelihood assigns 
the bulk of the fluctuations to CMB signal: Cl^cmb = 

72tl /iK, Co'lL < 20 /iK. As shown 
in table 9, these values are rather insensitive to the exact choice of the spectral index n. 

The two-component joint likelihood analysis on the IF data places the bulk of the 
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signal on the atmospheric component: CQatm — 189^39 and Cq(.^^ = 159l!lg| /xK, which 
in flat band power estimate becomes l)C^/(27r) = 4.1+^1 • 10"^. Treating the 

cahbration uncertainty as a systematic effect, the above results for the CMB component 
become : Cl[^^^ = 159lp//K and yje{e + 1)Ci/{2tt) = 4.lt|t • 10"^ Having a large 
atmospheric signal in our IF data is expected because the IF demodulation is known 
to be less efficient in removing linear gradients in the atmospheric signal. However, a 
significant fraction of the signal is projected to the CMB component which is marginally 
consistent with the 2F result. When testing for the presence of Galactic contamination the 
two-component likelihood analyses yield low significance detection for a dust component, 
and only upper limits to contamination by free-free or synchrotron emission. 



8. CONCLUSIONS 

This work demonstrates that it is possible to achieve sensitivities of a few tens of /iK 
in the study of CMB temperature anisotropy at millimetric wavelengths from ground based 
observatories. This requires sites with stable and dry atmosphere and the combination of 
long observing periods, measurements at several frequencies and a careful subtraction of 
the atmospheric contribution. Other conclusions of this work are: 

1. The sensitivity achieved in the 2F demodulation allows us to identify a common signal 
between our channels with a value of = 72I24 A*K corresponding to a band power 
estimate of l)Q-/(27r) = 2.0l J;^ • IQ-^ at Z = 531?^ at 68% CL including the 
systematic effect due to calibration uncertainty. We believe calibration uncertainty 
should be treated as a systematic error and not added in quadrature to the error bar. 

2. Our value is consistent with the detection reported by the Saskatoon experiment at 
the same angular scale. 
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3. Two alternative arguments (one based on extrapolations from Galactic templates, 
the other on results from other experiments) allow us to discard a Galactic origin of 
the detected signal. An extension of the likelihood analysis allowing the presence of 
a signal with the assumed spectrum for the atmospheric emission also discards the 
possibility of the signal being caused by correlated atmospheric residuals. 

4. We have used a simple technique for reducing the atmospheric noise in millimetric 
observations. We have presented tests on the performance of this technique by looking 
at the power spectra of the data before and after applying it. The reduction of the 
noise is also observed by comparing the distribution of the noise before and after 
cleaning the data. 

5. We have detected the Galactic Plane crossing in all cleaned channels at both 
demodulations. This constitutes the best check of the performance of the technique. 

6. In all channels the IF data show an excess of signal with respect to the 2F data. This 
excess is suspected to be of atmospheric origin. The two-component joint likelihood 
analysis assigns a portion of the detected signal to CMB fluctuations at a level of 
Cq/^ = or in band power estimate ^/i{l+l)cJ{2^ = 4.ll|^ • IQ-^ at 
68% CL including calibration uncertainty as a systematic effect. Our detected signal 
in the IF demodulation corresponds to an £-range never before observed. 

We expect to improve the results presented in this work with the addition of the new 
data taken during the Spring of 1996. 
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Fig. 1. — Window function for the IF and 2F demodulations at angular separation ■0=0. 

Fig. 2. — Astronomical calibration by observing the Moon at the IF demodulation. Solid 
lines are the data and dashed lines the predictions using the Lunar models cited in text. 

Fig. 3. — Meteorological conditions during the 1994 observing campaign. 

Fig. 4. — Power spectra of the in- and out-phase components at both demodulations and for 
all channels for a typical night of observation. Temperatures in this figure refer to antenna 
temperature values. 

Fig. 5. — Distribution of rms values for the whole campaign before applying the cleaning 
technique for the IF (solid line) and 2F (dash-dot line) demodulations. Note the height of 
the last bin in each histogram, as it contains the contributions of that bin plus all following 
bins. Temperatures in this figure refer to thermodynamic temperature values. 

Fig. 6. — Mean auto-correlation curves for all channels and both demodulations of the raw 
data. Dotted lines represent the ±lcr limits. 

Fig. 7. — The different stages of the cleaning technique. Panels in the left column show the 
raw data. Center panels display the output of equation (3), superimposed the sinusoidal fit 
(thick dashed line). Right panels show the cleaned data with the baseline removed. All plots 
have been brought to a bin size of 3° in RA for display purposes. Temperatures in this figure 
refer to thermodynamic temperature values. 

Fig. 8. — Final data sets in the regions to be considered for the posterior Likehhood analysis. 
The data have been binned to 3° in RA for clarity. We show the IF and 2F profiles indicating 
the instrument response to point sources. Temperatures in this figure refer to thermodynamic 
temperature values. 
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Fig. 9. — Power spectra of a typical night of observation in thermodynamic temperature at 
10s before (thin hne) and after (bold line) applying the cleaning technique. The dashed lines 
represent the upper limits to the instrument noise as estimated in section 4.1. Temperatures 
are expressed in thermodynamic units. 

Fig. 10. — Distribution of rms values after cleaning and editing for all the observing nights 
eligible to be cleaned. Notice the height of the last bin in each histogram for it contains 
the contributions of that bin plus all following bins. Temperatures in this figure refer to 
thermodynamic temperature values. 

Fig. 11. — Comparison of the recovered Galactic Plane crossings at channels 1, 2 and 3 and 
at both demodulations with the predictions in thermodynamic temperature. 

Fig. 12. — Mean cross-correlation curves between overlapping sections of individual nights 
used to generate the final data sets at any two channels. 

Fig. 13. — Contour plots of the hkelihood surface of the joint analysis on channels 1,2 and 
3 for the IF demodulation when a second component other than CMB is allowed. The four 
indeces represent the four relevant foregrounds: dust, atmosphere,free-free and synchrotron 
emission. Contour levels represent the confidence levels at 68% (solid hne) and 95% (dashed 
hne). The X symbol indicates the position of the hkelihood surface peak. Temperatures on 
the axis are in thermodynamic units. 

Fig. 14. — Contour plots of the likelihood surface of the joint analysis on channels 1,2 and 
3 for the IF demodulation when a second component other than CMB is allowed. The four 
indeces represent the four relevant foregrounds: dust, atmosphere,free-free and synchrotron 
emission. Contour levels represent the confidence levels at 68% (solid hne) and 95% (dashed 
line). The X symbol indicates the position of the likelihood surface peak. Temperatures on 
the axis are in thermodynamic units. 
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Table 1. Results From Fits To Extended Moon Transits. 



Chi 0.94 ±0.10 2.89 ±0.09 1.01 ± 0.27 2.80 ±0.29 

Ch 2 0.88 ±0.07 2.90 ±0.06 0.92 ±0.19 2.85 ±0.19 

Ch 3 0.85 ±0.06 2.90 ±0.05 0.88 ±0.17 2.86 ±0.16 

Ch 4 0.82 ±0.07 2.91 ± 0.06 0.84 ±0.13 2.88 ±0.13 
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Table 2. Correlation Between Channels Before Applying Atmospheric Correction. 



IF DEMODULATION 2F DEMODULATION 



Whole Campaign Final Data Set Whole Campaign Final Data Set 



Ch 1 - Ch 4 0.965 ± 0.010 0.990 ± 0.010 0.884 ± 0.023 0.91 ± 0.09 

Ch 2 - Ch 4 0.983 ± 0.008 0.995 ± 0.011 0.977 ± 0.007 0.991 ± 0.010 

Ch 3 - Ch 4 0.99963 ± 0.00010 0.99984 ± 0.00010 0.9956 ± 0.0011 0.995 ± 0.010 
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Table 3. Data Used To Generate Final Data Sets. 



IF DEMODULATION 




PtTCditagc 


NuuiIkt nights 


Mean rins (uiK) 


Amplitude Baseline (iiiK) 


Ch 1 


17.2 % 


16 


0.392 ± 0.024 


0.78 ± 0.19 


Ch 2 


13.3 % 


14 


0.48 ± 0.07 


1.23 ± 0.23 


Ch 3 


16.8 % 


12 


1.30 ± 0.14 


1.81 ± 0.23 


2F DEMODULATION 


Ch 1 


30.1 % 


24 


0.275 ± 0.024 


0.46 ± 0.08 


Ch 2 


15.4 % 


18 


0.233 ± 0.018 


0.85 ± 0.23 


Ch 3 


14.5 % 


15 


0.86 ± 0.12 


1.28 ± 0.25 



Note. — 



Mean rms and baseUne ampHtudes are expressed in thermodynamic units. 
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Table 4. Noise Spectrum {mK s^/^) Before And After Applying Atmospheric 

Subtraction. 





IF Out 


IF Out 


IF In 


IF In 


2F Out 


2F Out 


2F In 


2F In 


Channel 


1 Hz 


0.001 Hz 


0.001 Hz 


0.001 Hz 


1 Hz 


0.001 Hz 


0.001 Hz 


0.001 Hz 










CLEANED 








CLEANED 


1 


4.8±0.4 


6.9±0.6 


52±5 


8.1±1.0 


3.6±0.3 


6.7±0.6 


8.8±0.7 


6.1±0.7 


2 


1.13±0.10 


1.87±0.17 


115±10 


3.9±0.5 


1.04±0.11 


3.5d=0.4 


12.3±1.1 


3.3±0.4 


3 


1.99±0.18 


3.6±0.3 


570±50 


20.0±2.3 


1.81±0.17 


16.0±1.5 


55±5 


12.7±1.5 


4 


1.20±0.11 


18.9±1.7 


750±60 




1.07±0.10 


27.6±2.6 


94±8 





Note. — 



Mean rms and baseline amplitudes are expressed in thermodynamic units. 
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Table 5. Values Of pi4 For The Models Under Consideration. 



IF DEMODULATION 2F DEMODULATION 



Model 


Reference 


P14, 


P24 


P34 


Pl4 


P24 


P34 


1 


1 


0.1167 


0.1515 


0.6458 


0.1567 


0.1533 


0.6452 


2 


1 


0.1249 


0.1533 


0.6453 


0.1750 


0.1551 


0.6446 


3 


3 


0.1292 


0.1786 


0.6724 


0.1458 


0.1787 


0.6720 


4 


3 


0.1291 


0.1720 


0.6624 


0.1505 


0.1720 


0.6618 


5 


4 


0.1233 


0.1528 


0.6454 


0.1725 


0.1546 


0.6448 


6 


5 


0.1229 


0.1859 


0.6778 


0.1273 


0.1851 


0.6775 


7 


6 


0.1147 


0.1555 


0.6511 


0.1397 


0.1556 


0.6506 


8 


7 


0.1137 


0.1780 


0.6739 


0.1178 


0.1773 


0.6736 


9 


8 


0.1195 


0.1788 


0.6723 


0.1269 


0.1779 


0.6719 


10 


9 


0.1386 


0.2229 


0.7063 


0.1379 


0.2224 


0.7061 


11 


10 


0.1552 


0.1654 


0.6438 


0.1678 


0.2105 


0.6427 


12 


11 


0.1299 


0.1748 


0.6649 


0.1464 


0.1749 


0.6644 



References. — (1) Boulanger et al. 1996; (2) Reach et al. 1995; (3) Wright et al. 1991; (4) Bersanelli et al. 1995; (5) 
Davies et aL 1996a;(6) Kogut et al. 1996a; (7) Kogut et al. 1996b; (8) de Bernardis et al. 1991; (9) Banday & Wolfendale 
1991; (10) Page et al. 1990; (11) Fischer et al. 1995. 
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Table 6. Expected rms Values Due To Diffuse Galactic Emission. Units Of iJ,K. 



IF DEMODULATION 2F DEMODULATION 



Model 


Reference 


Description 


rmschi 


rmsch2 


rmschs 


rraschi 


rmsch2 


rmschs 


1 


1 


Z^" (XJ/2b„(17.5) 


1.2 


2.4 


12.3 


0.6 


1.3 


6.7 


2 


1 


oc 1.25^18.2) 


1. 1 


2.2 


11.3 


0.6 


1.2 


6.1 


3 


2 


oc 1.2(5^(20.4)+ 
+ 6.7 X B^(4.8)) 


1.5 


3.1 


13.6 


0.8 


1.7 


7.4 


4 


3 




1.5 


2.9 


13.1 


0.8 


1.6 


7.1 


5 


4 


oc i/2b^(18.0) 


1.1 


2.3 


11.6 


0.6 


1.2 


6.3 


6 


5 


oc i/1-4b^(23.3) 


2.5 


4.9 


19.6 


1.3 


2.6 


10.6 


7 


6 


1^ oc i/1-9B,(18) 


1.3 


2.8 


13.7 


0.7 


1.5 


7.4 


8 


7 


7^' oci/1-5Bi,(20) 


2.5 


5.2 


21.4 


1.3 


2.8 


11.6 


9 


8 


7^ oc i.1-5Bi,(22) 


2.2 


4.3 


18.2 


1.1 


2.3 


9.8 


10 


9 


Z^' oc iyB„{22.1) 


7.0 


12.5 


41.0 


3.8 


6.8 


22.2 


11 


10 


oc i/2B„(22.1) 


0.9 


1.5 


8.0 


0.5 


0.8 


4.3 


12 


11 


oci/1-6B,(24) 


1.6 


3.1 


13.5 


0.8 


1.6 


7.3 



References. — (1) Boulanger et al. 1996; (2) Reach et al. 1995; (3) Wright et al. 1991; (4) Bersanelli et al. 1995; (5) 
Davies et al. 1996a;(6) Kogut et al. 1996a; (7) Kogut et al. 1996b; (8) de Bernardis et al. 1991; (9) Banday & Wolfendale 
1991; (10) Page et al. 1990; (11) Fischer et al. 1995. 

Note. — 1^ and Bi, (T) stand for the dust spectrum and a black-body spectrum at a temperature of T Kelvin, respectively. 
Models (3) and (4) differ in the dust spectrum at low |6|. 
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Table 7. Basic Statistic Figures Of Final Data Sets. 



RAi RA2 





Range 


Mean 


Mean 


Mean 


Range 


Mean 


Mean 


Mean 




(°) 


# points 




rras {iiK) 


(°) 


# points 




rms i/J-K) 


Ch 1 IF 


[224,285] 


7.9 


130.6 


134.3 


[331,369] 


10.8 


103.1 


123.7 


Ch 2 IF 


[230,285] 


8.3 


118.9 


157.9 


[331,361] 


7.4 


113.2 


110.6 


Ch 3 IF 


[224,285] 


8.8 


312.6 


340.5 


[331,367] 


8.7 


310.6 


443.6 


Ch 1 2F 


[206,285] 


12.6 


74.7 


62.5 


[331,381] 


13.2 


74.1 


74.4 


Ch 2 2F 


[229,285] 


9.1 


68.5 


71.0 


[331,360] 


9.5 


53.3 


58.5 


Ch 3 2F 


[236,285] 


8.3 


236.4 


197.9 


[331,362] 


9.5 


187.7 


204.7 



Note. — Mean cr refers to the mean error bar associated with each in the final data set. Mean rms is the weighted rms 
along the indicated RA range. Both quantities expressed in thermodynamic units. 
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Table 8. Likelihood Results On Individual Data Sets. Values In /iK CMB. 



IF DEMODULATION 



2F DEMODULATION 



Channel 



(Co)V2 
RAi 



{Coy/' 
RA2 



{Co)'/' 
RA1+RA2 



{Co)'/' 
RAi 



{Co)'/' 
RA2 



{Co)'/' 
RA1+RA2 



Ch 1 
Ch 2 
Ch 3 



130 



+63 
-52 



7(34+219 



124 



+71 
-58 



155 



+83 
-58 



127 



+44 
-37 



155 

711; 



+83 

-58 



-168 
-135 



< 123 



91 



+47 

-43 



9oi:: + 158 
"3-181 



99 



+53 
-45 



132 



+77 

-56 



c;qi +284 



71 



+34 
-37 



106 



+37 

-33 



070+169 

"^'"^-163 



Note. — 



The stated confidence intervals do not include calibration uncertainties. 
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Table 9. Joint Likelihood Results. Values In CMB. 



IF DEMODULATION 2F DEMODULATION 



Channels 


(Co)V2 (^K) 


Cross-Corr i-j 


(Co)V2 (^K) 


Cross-Corr i-j 


Ch 1 X Ch 2 
Ch 1 X Ch 3 
Ch 2 X Ch 3 
Ch 1 X Ch 2 X Ch 3 


1291^^ 


0.19 ±0.07 
0.085 ±0.013 
0.31 ± 0.06 


79+35 
'^-37 

79+26 


0.149 ±0.026 
0.049 ±0.010 
0.26 ±0.05 



Note. — The stated confidence intervals do not include calibration uncertainties. 
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Table 10. Likelihood Results Allowing A Foreground Component. Values In /^K CMB. 



IF DEMODULATION 2F DEMODULATION 



n 


Foreground 


<l6 il^K) 








+2.0 


DUST 


128+41 






<6 


+1.5 


DUST 


128tg 




79+26 
'-^-22 


< 8 


+0.4 


ATMOSPHERE 


147t^^ 


150tf, 


72tlt 


< 16 


+0.2 


ATMOSPHERE 






79+26 
'-^-24 


< 18 


+0.0 


ATMOSPHERE 






79+26 
'-^-24 


< 20 


-1.8 


FREE-FREE 


U741 


< 177 


72tlt 


< 62 


-2.1 


FREE-FREE 




< 189 


79+26 
'^-26 


< 74 


-2.4 


FREE-FREE 


153^39 


< 201 


79+26 
'^-26 


< 86 


-2.7 


SYNCHROTRON 


153^39 


< 210 


79+26 
'^-26 


< 96 


-3.0 


SYNCHROTRON 




< 216 


79+26 
' ^-28 


< 106 


-3.3 


SYNCHROTRON 


156t^^ 


< 222 


79+26 
'^-28 


< 116 



Note. — The stated confidence intervals do not include calibration uncertainties. 
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